Three-body problem at finite temperature and density 



A. N. Kvinikhidze^' ^ and B. Blankleider^ 

^A. Razmadze Mathematical Institute, 
Georgian Academy of Sciences, Aleksidze Str.l, Tbilisi 0193, Georgia 
^Department of Physics, Flinders University, Bedford Park, SA 5042, Australia 

(Dated: February 9, 2008) 

Abstract 

We derive practical three-body equations for the equal-time three-body Green function in matter. 
Our equations describe both bosons and fermions at finite density and temperature, and take into 
account all possible two-body sub-processes allowed by the underlying Hamiltonian. 

PACS numbers: 21.45. +v, 21.65. +f, 24.10. Cn, ll.lO.Wx, ll.80.Jy 



1 



I. INTRODUCTION 



Three-body correlations play an important role in describing many aspects of many- 
body systems. In early studies of nuclear matter, where the main focus of interest was 
its binding energy, three-body correlations were found to contribute a small but significant 
part [H, 12^ 0, IJ More recently, it has been the study of in-matter three-body systems 
themselves, that has been of main interest. Indeed, the in-matter three-body problem plays 
an important role in describing a large variety of interesting phenomena in many-body 
systems. For example, in order to understand the formation of bound states in heavy ion 
collisions, three-body calculations are needed to study the modification of the binding energy 
and wave function of a three-nucleon bound state due to nuclear matter of finite density 
and temperature 0|. Similarly, studies of the binding energy of three quarks are of 
relevance to the understanding of color superconductivity and phase transitions in quark 
matter 0, |^. Three-body calculations are also needed to describe nonequilibrium processes 
of cluster formation in an interacting many-body system pl^Y ^ fundamental role 

in determining the two-particle-one-hole (pph) and two- hole-one-particle (hhp) contributions 
to the self-energy of the single-particle propagator 

The goal of the present paper is to formulate three-dimensional equations for the finite 
temperature in-matter three-body problem, that take into account all possible two-body 
sub-processes allowed by the underlying Hamiltonian. To put this goal into context, it is 
worthwhile to briefly review the progress made so far on this subject. From the very be- 
ginning, it was recognized that Faddeev's approach [l^ provided a powerful tool in the 
description of few-body properties in Quantum Mechanics. It is therefore natural that not 
long after its formulation, this approach was also applied to quantum field theory, first within 
a four-dimensional formulation [13|, and then within a three-dimensional one obtained by 
equating times in four- dimensional Green functions These early formulations were for 
three particles in vacuum. With the application of quantum field theoretical methods to 
statistical physics 1^ 3], it became possible to apply Fadeev's approach also to the field- 
theoretic description of three particles within a many-body environment. However, one ma- 
jor obstacle in formulating a practical description in this way, is the hole contribution to the 
single-particle propagator in the form of an advanced part [see, for example, Eq. (j^], which 
is not present in the quantum mechanical description of particles in vacuum. The presence 
of this hole contribution makes the field-theoretic description inherently four dimensional, 
even in the non-relativistic case. The first steps in applying the Faddeev approach to the 
many-body enivornment avoided this problem by utilizing a Bethe-Goldstone type of mod- 
ification of the Faddeev equations, which involves of a simple momentum cut-off restricting 
the intermediate state particles to be above the Fermi surface Although such modified 
equations can be treated with the Faddeev method, they do not take into account the hole 
contributions which reside in the advanced parts of the single particle propagators. The 
way beyond this approximation was proposed by Schuck, Villars, and Ring Ij], who used 
equal-time Green functions in order to obtain a three-dimensional field theoretic description. 
To derive their equation for the zero-temperature equal time three-body wave function, they 
approximated the effective pair-interaction kernels by terms linear in the physical two-body 
potentials. Since the exact expression for the effective pair-interaction kernel involves an 
infinite series of higher order terms as well [see Eq. ()48p]. the linear approximation cannot 
be considered as satisfactory for the strong coupling case, e.g., when two-body bound states 
are possible. The current state-of-the art formulation [islj, which has been used extensively 
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for calculations jE B 0, 0, S 113, 113 be considered as the model of Ref. 17 1 extended to 
finite temperatures, with the extension being performed using the imaginary time formalism 
of perturbation theory 0|. 

In this context, the goal of the present paper is to formulate practical field-theoretic three- 
body equations, valid at finite temperature and density, that take into account the whole of 
the above mentioned series for the effective pair interaction. Like Refs. fl^, we use 

equal-time Green functions to formulate a three-dimensional field theoretic description. We 
show that Faddeev's idea which renders the three-body kernel compact, namely, to reexpress 
the three-body equations in terms of two-body t matrices rather than two-body potentials, 
also enables one to sum up, exactly, the infinite series of Eq. ()48|) for the pair-interaction 
kernel. 



II. IN-MATTER FOUR-DIMENSIONAL THREE-BODY EQUATIONS 

The interactions of three identical particles at finite density and temperature are described 
in quantum field theory by the Green function Q defined by 

{27iy6^{p[ + P2 + ^3 - Pi - P2 - P3)Q{p'lP2P3' P1P2P3) = j rf^?/lrf^?/2'^^l/3C?''a;iC?^X2«i^X3 

eiK-yiM-s/2+p^-2/3-Pi xi-p2.x2-p3-x3) [pT[m{yi)m{y2)^{y3)^\xi)^!\x2)^\xi)]] (1) 

where \l/ and are Heisenberg fields with respect to the Hamiltonian K = H — fiN, T is 
the time ordering operator and 

P=T^ (2) 

is the statistical operator of the grand canonical ensemble Besides being the central 

quantity for the description of three-body observables, this Green function is also needed to 
calculate the vacuum properties of the system with the help of the dressed single particle 
propagator; for example, in the four-point interaction model, the single particle self-energy 
diagram is completely defined by the particle-particle-hole (pph) Green function jlQ^ . 

In the zero temperature case {(3 = l/fc^T oo) only the ground state survives in the 
trace of Eq. (0), and the Green function reduces to the usual QFT expectation value 

{2tcY6^{p[ + p'2 + p'3 -P1-P2- P3)Q{p'iP2P3'^ P1P2P3) = j d'^yid'^y2d^y3d^xid^X2d^X3 

^M-yi+p'ry2+p'ry3-Pi-xi-P2-x2-P3-x3) (0|r^(|/i)^(|/2)^(|/3)^^(x3)^^(x2)^^(xi) |0) (3) 

where |0) is the physical ground state. For the latter, straightforward use of Wick's theorem 
gives a perturbation theory with Feynman rules. Analogously, two types of perturbation 
theory, so-called "imaginary-time" and "real-time", have been derived for Eq. [2o| . 
Here we shall use the real-time formulation of perturbation theory in which the number of 



degrees of freedom is doubled [2l|; this complication, with respect to the zero-temperature 



case, comes from the sum over the complete set of states (trace) in Eq. ((TJ).^ For example. 



^ A similar discussion on the basis of the imaginary-time formalism and its comparison with this note will 
be presented elsewhere. 
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the free one-body propagator, in the nonrelativistic case, is given by [2 



Hp) ^ Hp) 



p" — Lu + ie — uj — ie 



(4) 



where u = Up = p^/2m — fi (we take ^ = 1), /i is the chemical potential, and n, n are 2x2 
matrices whose elements are simple functions of the distribution function 

/(-) = 

with the upper sign (+) for fermions and the lower sign (— ) for bosons (see Appendix A). 
Correspondingly, the elementary vertices have an extra double- valued index for each particle 
leg; e.g., the four-point interaction v used to define the Hamiltonian of Eq. (fTTj) . enters the 
formalism with doubled degrees of freedom through the quantity v whose matrix structure 
is given as Vij^i = (^ii^ji'^fci'^n — Si2Sj2Sk2Si2)v- Note that the first diagonal element is just 
the potential itself, vim = v. Similarly, the first diagonal element of Eq. d(i, for which 

nn(p) = 1 - nn(p) = ±;^^ = ±/(^). (6) 

corresponds to the usual one-body free Green function at zero temperature In some 
cases one can use a simplified propagator consisting of just this dd element of Eq. (jU [23l |. 

For the identical particle case considered here, the field theoretic expression of Eq. ((T} 
automatically guarantees the appropriate symmetry of the three-particle Green function 
Q. Moreover, in the doubled degrees of freedom formalism, the matrix Green function G 
whose first diagonal element is ^, is likewise properly symmetric in the case of bosons, and 
antisymmetric in the case of fermions. On the other hand, the disconnected Green function 
Go defined by 

Goip'M^PM) = d{p,) d{p,) d{ps) {2t:)'5\p', - P2) {2t:)H\p', - ps), (7) 

where d{pi) is the dressed propagator of particle i, does not possess identical particle sym- 
metry; thus Go is not equal to the fully disconnected part of G (which we shall denote 
by Grf). Indeed, it can be easily shown that to obtain G^, one need only symmetrize (or 
antisymmetrize) Gq according to the equation 

Go(l'2'3', 123) = G^(l'2'3', 123) = Gd(l'2'3', 123) (8) 

p 

where the sum is over all permutations P of either the initial or final state particle labels, 
and for fermions, is understood to include a factor (—1)''^ = +1 or —1 depending on whether 
the permutation is even or odd, respectively. In Eq. (jHJ we use a symbolic notation where 
integers represent the momenta plus all quantum numbers of the corresponding particles, 
with primes distinguishing the final states. 

Defining the kernel K to be the set of all possible three-particle irreducible Feynman 
diagrams for the 3 — 3 process, we may write the equation for the Green function G as ^2j| 

G = G^ + jfoKG (9) 
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where the 1/3! factor reflects the fact that both G and K are fully symmetric or antisym- 
metric in their particle labels. The disconnected part of K, indicated by subscript d, can be 
expressed in terms of the identical particle two-body potential v as 

Kd{l'2'3', 123) = ^(2'3'> 23)d'\l)6{l' , 1) (10) 

LcRc 

where 5(1', 1) represents the momentum conserving Dirac 6 function {2tt)^6^{p'i —pi), while 
Lc and Rc indicate that sums are taken over cyclic permutations of the left labels (1'2'3') 
and right labels (123), respectively (note that the sums are restricted to cyclic permutations 
because the potential v is already properly symmetric or antisymmetric in its labels). 
Defining 

V,{1'2'3,123) = v{fk',jk)d~\z)6{i',i) (11) 
where {ijk) is a cyclic permutation of (123), we have that 

i^d = 5^(^l + ^2 + V^3), (12) 

Pc 

where it makes no difference over which labels, left or right, the cyclic permutations are 
taken. Denoting the connected part of the kernel by Kc, we define the 3—^3 potential V 
by 

^=^(^l + ^2 + ^3) + ^^c. (13) 

2 

Although V is not fully symmetric or antisymmetric, it does have the useful symmetry 
property 

P,,VP., = V (14) 

where Pij is the operator that exchanges the i'th and j'th momentum, spin, and isospin 
labels. Since 

K = J2V, (15) 
p 

Eq. Q can be written as 

G = + GoVG. (16) 

Formally, Eq. ()16p differs from the equivalent relation for distinguishable particles, G = 
Go + GqVG, only in the explicit symmetrization of the inhomogeneous term. 

For the sake of simplicity, we consider the model given by the second quantized Hamil- 
tonian 

H = 'Y^uia\ai + 'Y^v{1234:)a\a2a4a3 (17) 

1 1234 

where Ui = q^/2m is the single-particle non-relativistic kinetic energy, and w(1234) = 
^(siqi, S2q2, Ssqs, S4q4) is a Galhlei invariant function of spins Si and three-dimensional 
momenta qj of particles 1, 2, 3 and 4. Then the first diagonal element of the pair interaction 
potential v of Eq. (fTTj) . fim, is not v in the general case of finite density and temperature, 
as the presence of the medium makes it possible for particle-hole pairs to be exchanged in 
the t-channel. As a result of these exchanges v = v{siqi, S2q2, Ssgs, 84^4) also depends on the 
energy variables, q^, g^, ^3, Qa (of which only three are independent since qi+q2 = Qs + Qa)- In 
the case of zero density and temperature, /x = 0, T = 0, we obtain pure quantum mechanics 
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with particle number conservation and fun = v. In this case, single particle propagators 
do not have an advanced part, n = 0, and this fact allows one to derive the standard 
three-dimensional Faddeev equations for the three-particle system where input two-body t 
matrices (in three-body space) are obtained from physical two-body t matrices (in two-body 
space) simply by subtracting the spectator energy from the total energy. In the next two 
sections we discuss how three-dimensional three-body equations can also be derived for the 
case of finite temperature and density. 



III. THREE-DIMENSIONAL EQUAL-TIME REDUCTION 

Being a four-dimensional integral equation, Eq. is not very convenient for practical 
calculations - it involves integrations over relative times (or relative energies) which, because 
of the presence of cuts and singularities in all four quadrants of the integration plane, 
cannot be easily handled numerically. For this reason, we would like to implement a three- 
dimensional reduction of this equation. To do this we follow the current literature and effect 
this reduction by equating times in initial states, and separately, in final states. A quantity 
A with such equated times will be denoted by i{A) (an extra factor of i has been included 
for later convenience). Thus, our central quantity is the two-time Green function {Q) which 
is obtained from the four-dimensional Green function Q by equating times as just described. 
In momentum space, {Q) is therefore given by 

{2Txf5^{Y>'i + P2 + Ps - Pi - P2 - Pa) ^{Q){E, p'lPaPs. P1P2P3) 

d^yi rf^2/2 d'^xi d^X2 ci^xs dt e^(^*-p'i yi-P2 y2-pJ,-y3+Pi-xi+P2-x2+P3-x3) 

Tr {pT [M/(t, y^)^it, y2)^{t, y3)^^(0, X3)vl/t(0, X2)vl/t(0, xi)] } . (18) 

In the doubled degrees formalism, the two-time Green function of Eq. (jl8|) is the first 
diagonal element of the matrix Green function (G). The goal of this section is to develop 
the three-dimensional integral equation, analogous to Eq. (fT^. for {G). Previously, such 
three-dimensional equations for the two-time Green function have been considered in the 
context of (zero density) relativistic quantum field theory in Ref. ji^ for the two-particle 
case, and in Ref. for the three-particle case. We shall base our derivation of the in- 
matter three-body equation upon that of Ref. 3|. By contrast, other recent many-body 



lany- 

formulations of three-body equations have been closely related to the work of Ref. [17[ . The 
difference between these two approaches will be discussed in due course. 

It can be shown that the momentum space two-time Green function of Eq. (jl8|) can be 
obtained directly from the four- dimensional Green function of Eq. (Q) by integrating out all 
the relative energies: 

■/r\fjp I ' ' \ [ dp'i dp'2 dp'^^ dpi dpi dpi , , , 

i{g) [E, P1P2P3, P1P2P3) = J ^ ^ 27 27 27 ^^^i^2P3' P1P2P3) 

(27r)2 5(p;° + + p' - E) 5{pI + + p° - E). (19) 

Correspondingly, in the doubled degrees of freedom formalism, the equal time matrix Green 
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function (G) is related to its four- dimensional counterpart G by 

• /r-v/TT. / / / A f dp'i^ dp'^^ dpi dpi dpi , 

{2nf 5{p\' + p',' + j9'3° - E) 5{p\ +pI+pI-E). (20) 
Thus the four-dimensional free three-body Green function, defined as 

Gi{p\p'2P',,ViV2Pz) = d^iPi) d^{p2) d^iPs) {27rf 6\p', - p^) 5\p', - ps), (21) 
where d^ is given by Eq. (@]), leads to the following two-time free three-body Green function: 

(G'^)(E,p;p'2p'3,PlP2P3) = (Go0(i5,PlP2P3)(27r)«53(p'2-p2)53(p'3-p3), (22) 

where^ 

(GoO(i5,PiP2P3) = -^ / ^^df{E-pl-plY>i)df{plv2)df{plY^,) 

^(Pl)^(P2)^(P3) ^ ^(Pl)^(P2)^(P3) ^23) 



E -ujp^- tUp^ -Up^+ie E - cup^ - Wp^ - uj^.^ - ie 

In a similar way, we can define the Green function T> of two identical particles at finite 
density and temperature as 

(27r)V(pi +P2 -Pi - P2)1^ip[P2,PiP2) = J d'^yid'^y2d'^Xid'^X2 

^^{p'ryl+P'2■V2-Pl■xl-P2■x2) |pr[^(?/i)^(2/2)^^(x2)^^(xi)]} (24) 
from which the two-time Green function (P) follows: 

{2nf 5{p\' + p',' - E) 5{p\ +pI~E). (25) 

In the doubled degrees of freedom formalism a similar equation relates the two-body matrix 
Green function D to its two-time version {D). Applying this relation to the free two-body 
Green function 

dUp'iP2,PiP2) = d^iPi) df{p2) {2nY6\p', - p^), (26) 
one obtains the two-time free two-body Green function: 

(Di()(E,p;p^,PiP2) = (DoO(^,PiP2) i2nf5%p', - P2), (27) 



where 



{D|;)iE,p^p2) = -^ J df {E - p^p,) d^ {plp^) 



'^lPij^lP2j n[pi)n[jp2] 



E — a;p^ — a;p2 + ie E — u^^ — uo^,^ — ie 



(28) 



^ To save on notation, we use the same symbol with differing numbers of arguments to represent disconnected 
Green functions with and without the momentum conserving 5 functions. 
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It also convenient to express the one-body propagator as 

{df)iE,p) = -^dfip)l,_^ = "^P^ . + "(P^ . . (29) 

--^ E — Up + le E — uj-p — le 

A. Three-body equal time quasi-potential 

One can write a three-dimensional equation for the two-time Green function (G) of the 
same form as four-dimensional Eq. (jlfij) : 

(G) = (Go)^ + (Go)\/(G), (30) 



where V is the three-dimensional "quasi-potential". As is well known j25|, V should be 
expressible in terms of a perturbation series in the four-dimensional potential V\ however, 
this series involves the inverse of the disconnected Green function (Go)) which in the many- 
body case is a problem, as this inverse may not exist. To see this, it is sufficient to consider 
the free propagator of Eq. for the case of fermions at T = where real time perturbation 
theory is expressed in single degrees of freedom. In this case the functions n and n reduce 
to 

n(p) = e{pF - |p|), n(p) = 1 - n(p) = ^^(IpI - pf) (31) 

where pp = fJ'{0) is the Fermi momentum. It is clear that n and n are projection operators, 
and this in turn means that (Gg) is also a projection operator. Thus (Gq) does not project 
onto the full space of momenta, but only onto the sub-space projected by the operator 

Af = n(pi)n(p2)n(p3) + n(pi)fi(p2)fi(p3); (32) 

as a result, (Gq) cannot be inverted in the full space of momenta. Exactly the same con- 
clusion is reached in the case of finite temperature in the formalism with doubled degrees of 
freedom, as n and n remain projection operators, but now as matrices. More specifically (see 
Appendix A), n and n satisfy the following relations which define their projection properties: 
n + n = g,ngn = ngn = 0, ngn = n, and ngn = n where g is defined as 

(7 = ^ Q ^ ^ for fermions, ^ ~ '^l ^ bosons. (33) 

To get around the non-invertability of (Go), we introduce a modified Green function 

Go = {Go) + iggg-Af)A (34) 

which for non-zero A is not singular and can be inverted. In Eq. (j34p . ggg is a direct product 
of g^s where one g is in the space of particle 1, another is in the space of particle 2, and the 
third is in the space of particle 3. It is_ important to note that Go is identical to (Gq) in the 
subspace projected by Af, i.e., MgggGo = Afggg{Go) = (Gq). The operator A is required 
to be fully disconnected, so that 

Go{E, p[p',pi„ P1P2P3) = Go{E, P1P2P3) {27rY6\p', - P2) 6'{p', - Pa), (35) 

but can otherwise be chosen according to one's own convenience (later we will also be able 
to let A go safely to zero). For the free case, one can write down the inverse of Gq explicitly: 

{GI)-\E, P1P2P3) = 999 [^{E - Up, - Up, -Up,) + {ggg - ^)A-^] ggg. (36) 
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As a consequence of introducing Gq, we shall likewise introduce a modified full three- 
dimensional Green function G defined as 

G = {G) + {ggg-Af)A'' (37) 

and redefine the three-dimensional quasi-potential V to satisfy the equation 

G = G^ + GqVG. (38) 

As it stands, Eq. ()16|) is not convenient to work with as the inhomogeneous term con- 
tains an explicit sum over the permutations of Go's labels. For this reason we define the 
unsymmetrized Green function G" as the solution of 

= Go + GoVG"". (39) 

The full Green function G is then obtained by summing over the permutations of the 
right-hand labels of G", which we symbolically write as G = G"P. Similarly, in the three- 
dimensional case we define 

G" = Go + GoV^G'' (40) 

where G = G^P. It follows that 

G" = (G") + (^?(7(7-Ar)A. (41) 

Iterating Eq. (j39p . equating initial and final times, and using Eqs. (j34j) and (j4H) . we obtain 

G" = Go + (Gol^Go) + (GoFGol^Go) + . . . (42) 

where in momentum space the angle brackets indicate the integration of relative energies as 
in Eq. ()2U|). Since the inverse Gq^ exists by construction. 

Go 'G" = 1 + Go^ {GoVGo) + Go ' (GoI^GqI^Go) + . . . (43) 

so that 

(G")-iGo = 1 - Go ^GoVGo) - Go '(Go^Go^Go) 

+ Go '(GorGo)Go '(GoV^Go) - . . . (44) 

Using this in Eq. (j4Up gives an explicit perturbation series for the quasi-potential: 



V = Go ^ 



(GoV^Go) + (Gol^Gol^Go) - (GoVGo)Go '(GoVGo) + . . . Go ^ (45) 



In the present case of a three-particle system, V consists of a sum of pair interactions and 
three-body forces, as given in Eq. (fT^. The quasi-potential V must be expressible similarly 

as 

V=^{V, + V2 + V^a) + ^k, (46) 

where is a pair interaction with particle i as spectator, and Kc is a three-particle irreducible 
connected term (the three-body force). We can thus write Vi as 

^i(p'iP2P3, P1P2P3) = (27r)353(p^ _ p,)l>,(p;.p'^, p.pfc; Pi) (47) 
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where, to save on notation, we have used the same symbol Vi, but with different arguments, 
to denote the pair interaction with and without the spectator 6 function. Note that, in 
general, the Vi without the 6 function still depends on the momentum pi of the spectator. 
It follows from Eq. (|45|) that 



(Co— Go) + (Cq— Gq— Go) — {Go-^Go)Go^ {Gq-^Go) + 



Go'- (48) 



We are interested in the usual case where the three-body force is neglected. Then Vi 
provides the sole interaction that describes three-body observables. However, what enters 
the three-body Faddeev equations is the pair-interaction t matrix in three-body space, Tj, 
defined in terms of the quasi-potential Vi by the Lippmann-Schwinger equation 

f. = V + ^VGof,. (49) 

In correspondence with Eq. (jlTj) . we shall write the disconnectedness structure of Tj as 

Ti(p'iP2P3, P1P2P3) = (27r)^5^(pi - Pi)Ti(p^.pfc, pjPk, Pi). (50) 

The task of constructing Tj appears to be formidable. Just to determine Vi by summing the 
infinite series of Eq. fl48|l seems already a practical impossibility. That is why most, if not 
all, works on this subject keep only the linear term in the input two-body interaction jl7l |. 
In our case, this would mean keeping only the first term of the series in Eq. fl48|) . However, 
as we shall show subsequently, there is another way of solving this problem which gives 
the exact t matrix Tj, namely the one that results from a complete summation of Eq. fl48|) 
followed by an exact solution of the Lippmann-Schwinger equation, Eq. ()49p. 

We finish this_subsection by pointing out that the pair interaction part of the three-body 
quasi-potential, Vj (without the spectator 6 function), is related in a most non-trivial way to 
the corresponding two-body quasi-potential Vi arising from equating times in the two-body 
Green function. Indeed if we write the two-body equivalent of Eq. fll6|) . 

D = Di; + ^DovD, (51) 

where D is the two-body Green function originating from Eq. ()24|) . Dq is the two-body 
disconnected Green function given by 

Do{p[P2,PiP2) = d{p,)d{p2){2'Kf5\p', - P2), (52) 

and V is the fully symmetric (or antisymmetric) two-body kernel, then the two-body quasi- 
potential V can be determined by following the same procedure as the one above to determine 
V. Clearly, we shall obtain the following perturbation series for v: 



2 

where 
with 



(Do-Do) + (Do-Do-Do) - {Do-Do)D,\Do-Do) + . . .J ^o ' (53) 
^o = (^o) + (^7^7-Ar(2))A(2) (54) 



Ar(2) = n(pi)n(p2) + n(pi)n(p2), (55) 

and A*^^^ is any convenient disconnected function that makes Do invertible. It is easy to 
see that in general, Vi will be related to Vi in a complicated and non-linear way through an 
infinite series. 
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B. Exact three-body equal-time disconnected t matrix 



As the determination of the exact quasi-potential Vi is difficult, as discussed above, it 
would seem that the exact determination of the corresponding t matrix Tj, needed as input 
to the three-dimensional Faddeev equations, is not attainable, at least not through the 
iteration of the exact Vi. Yet it often happens that working with t matrices directly is much 
simpler than working with the underlying potential - one example being the problem of 
disconnectedness in the three-body problem, which Faddeev solved by formulating equations 
in terms of two-body t matrices rather than two-body potentials. 

Thus, rather than expressing Vi in terms of Vi, as we have done above in Eq. ()48|) . we shall 
now attempt to express the exact Tj directly in terms of the corresponding four- dimensional 
disconnected t matrix Tj. We start by noting that even though Vi may not be known, one 
can nevertheless formally write down the exact as the solution of the Lippmann-Schwinger 
equation, Eq. (jiUI) . Thus, if is the disconnected part of with particle i as spectator, 
so that, by Eq. (HHll . 

Gr = Go + ^GoVG'l, (56) 

it follows that ^ 

= Co + -^GoTiGQ. (57) 

Defining the equal time three-body t matrix T through 

G" = Go + GqTGq, (58) 

it's clear that one can write 

f = ^(f,+f2 + n^+t (59) 

where the part of T which is fully connected. 

In the four-dimensional sector one can similarly express the unsymmetrized Green func- 
tion G" in terms of the three-body t matrix T as 

G" = Go + GoTGo (60) 

where 

T = V + VGoT. (61) 
By analogy with Eq. |TT|) and Eq. we write 

T=^{T, + T2 + T;)+T, (62) 

where 

Ti(l'2'3',123) = t{j'k',jk)d-\i)6{i',i) (63) 

is the disconnnected part of T with particle i as spectator and t{j'k',jk) = U is the four- 
dimensional two-body t matrix, and Tc is the connected part. It follows that 

T, = V + ^ViGoT,. (64) 
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The disconnected part of G", with particle i as spectator, is given by 



Gl = Go + \GoT,Go. (65) 



Performing the equal time operation on Eq. (j65|) . 



(G'n = (Go) + ^(GoT,Go), (66) 
it follows from Eq. ^ and Eq. (gH) that 

Gr = Go + ^{GoT,Go). (67) 

Comparing this equation with Eq. (j57|) . we obtain 

GoT,Go = (GoT.Go). (68) 

It is seen that in contrast to the quasi-potential Vi which is related to the four- dimensional po- 
tential Vi in a very complicated way, the t matrix Tj corresponding to the quasi-potentional, 
defined by the exact solution of Eq. (jUJ, is connected to the four- dimensional t matrix Tj 
in a very simple way. 

Indeed, using Eq. and Eq. (jUHj) . one can write the RHS of Eq. (jUHj) as 

(GoT.Go) = -^/^^^^^^ [Do.UDoW,PlP.P^)i27rr5\p[ - p.)dM) 
{2nf 5{p\' + p'/ + V'' - E) 5{p\ +pI+pI-E) 

i2nf 5{p.' + p'; + p',' - E) 5{p^ + p] + pI - E) 

dp^ 



{27rf6'{p[ - p,) I '^t{D,,UD,,){E - pI p;.p',, p.-p^) (rf,)(p°, p,) (69) 



where 



[DoiUDo.^ {p'jp'^, PjPk) = d{p'j)d{p'f^) t{p'jp'f^,pjPk) d{pj)d{pk) . (70) 
The result of Eq. can thus be written as 

GofiGo = {2-Kf5\p[ - Pi) {D^ADoi) ® (rf,), (71) 
or, without the spectator 5 function, as 

GoTiGo = {D^^UD^,) ® {di) (72) 
where the symbol ® denotes the convolution integral: 

i f°° 

a ® b{E) = — a{E- z)h{z) dz. (73) 
2vr i_oo 
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The above analysis of three-body Green functions can be repeated for two-body Green 
functions, thereby yielding the two-body version of Eq. ()68p: 

DoiDo= {DotDo), (74) 

where 

i = v + -vDoi. (75) 

Using this in Eq. ()72|1 we obtain the essential result 

Gof,Go = DoiUDoi®{di) (76) 

which expresses the exact t matrix Tj, forming the 'spectator plus interacting pair' input 
to the three-body Faddeev equations, in terms of a convolution of the spectator propagator 
and the sub-system equal-time two-body t matrix. 



IV. IN-MATTER THREE-DIMENSIONAL THREE-BODY EQUATIONS 
A. General description 

The unsymmetrized three-body equal time Green function G^ is given in terms of the 
three-body t jnatrix T by Eq. (|KHjl . If the three-body force Kc of Eq. pHjl is neglected, one 
can express T in the Faddeev form 



T = ^X, (77) 



i=l 



where 



X, = ^f, + ^f,,J2GoXk. (78) 



2 2 



Alternatively, one can express in terms of Alt-Grassberger-Sandhas (AGS) amplitudes 

= GAj + G,Ui,G, (79) 
where the Uij satisfy the AGS equations 

Uij = Gq ^6ij + - SikTkGoUkj- (80) 

k 

In either case, the input consists of the disconnected amplitudes Tj which are specified in 
terms of equal-time two-body t matrices U according to Eq. (fTBj) . 

The above equations constitute our general formulation of the finite temperature equal- 
time in-matter three-body problem. What is noteworthy is that the neglect of the three- 
body force Kc is the only approximation made; in particular, the four-dimensional two- 
body potential Vi of Eq. (fTTj) . which is specified by the underlying Hamiltonian, is included 
exactly and to all orders within the equal-time approach. Although the above equations 
can be used directly for calculations, in the next subsection we shall show that for the 
case of instantaneous potentials and effective single-particle dressings, they can be greatly 
simplified. 
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B. Description for instantaneous potentials and free-like dressed propagators 



Here we consider the commonly used approximations where the four-dimensional two- 
body potential v is assumed to be instantaneous, and where single particle dressings are 
taken into account only through effective masses and effective chemical potentials. 

Thus we assume that the dressed propagator d has exactly the same structure as the free 
propagator (i-^, given in Eq. but with a modified (effective) mass m* and a modified 
chemical potential jj,*. From Eq. ()29p this means we can write for particle i, 

{d,) = n,dl + n4t (81) 

= + «) (82) 

where 

dT = , (83) 

d1 being specified with +ie and with — ze, fii = n{pi), rii = n{pi), and uji = pf/2m* — /i*. 
The disconnected equal-time two-body propagator, (-Do); will then take the same form as 
(Dg) given in Eq. ()28|) . and we similarly write 

{Doi) = UjUkDl - n.nkDI (84) 
= (^S.) + (^0.) (85) 

where 

DT = . (86) 

E — ujj — tUk ie 

The disconnected equal-time three-body propagator, (Go)? can likewise be written as 

(Go) = UinjUkG^ + niUjUkG"- (87) 
= (GS) + (GS) (88) 

where 

G'-''^ = . (89) 

E — Ui — ujj — Uk ± le 

At the same time, the assumption of an instantaneous two-body potential means, in 
momentum space, that potential f , and therefore the corresponding t matrix t, do not 
depend of zero components of relative momenta. Thus {DotDo) = {DQ)it{DQ), so that 
Eq. becomes 

DoiDo = {Do)it{Do). (90) 

Now, because the disconnected equal-time two-body propagator (Dq) has the form specified 
by Eq. it follows that (Dq) = {Do)ggAf^'^\ and therefore (Do) = DoggAf^'^l Applying 
this to Eq. (jHOI), one obtains the useful identity 

{Do)i{Do) = {Do)it{Do) = DoiDo. (91) 

Interestingly, the quasi-potential v and the instantaneous potential v obey a similar equation, 

{Do)d{Do) = {Do)tv{Do) = DovDo. (92) 
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Indeed, for the case where v is instantaneous and {Dq) is specified by Eq. (IS^ . Eq. (fTI)) 
holds also for potentials, i.e., 

b^vb^ = (DovDo). (93) 

To see this, consider the second order term in v of Eq. under the assumption of instan- 
taneous V. Up to a constant factor, we have that 

{DqvDovDo) - {DovDo)b^'{DovDo) 

= {Do)iv{Do)tv{Do) - {Do)tv{Do)b,\Do)iv{Do) 

= {Do)tv{Do) (l - b,\Do)) iv{D^). (94) 

Then taking {Dq) to be given as in Eq. (jini), so that (Dq) = {Do)ggAf^^^ = boggAf^^\ 

{DovDqvDq) - {DovDo)b,\DovDo) 

= {Do)iv{Do) (^ggAf^^^ - b,'{Do)) iv{D,) (95) 

= {Do)tv{Do)b,' {boggM^^^ - (Do)) tv{Do) = 0. (96) 

In a similar way, all higher order contributions in v are zero in Eq. (j53p . and the result of 
Eq. follows. 

Applying Eq. (|^ to Eq. (fTUj) . one obtains 

Gof,Go = {Do,)U{Do,)(»{d,) (97) 

which is to be used as the input to the equal-time Faddeev equations discussed above. 

1. Split form of equations 

Under the assumptions of this subsection, the disconnected two- and three-particle 
dressed propagators, (Dq) and (Go); have the same projection properties as the correspond- 
ing free propagators. These projection properties lead to a substantial simplification of 
the in-matter three-body equations. We shall show that for instantaneous potentials and 
the free-like dressed propagators of Eq. (jHUj) and Eq. (jH^ . the in-matter Faddeev equation, 
Eq. ffTSj) . can be split into two equations, one involving only retarded parts of propagators, 
the other involving only the advanced parts. 

In view of Eq. (|97p . we begin by defining the quantities Tf^ and T^^ by^ 

GoT;^Go = (27r)353(p/ _ (Do.)t.(Do.) ® «), (98a) 

Gof^^Go = (27r)353(p/ _ p^) {DoMDo^) ® (rf^, (98b) 

so that 

f, = f^ + f^. (99) 



^ The superscripts of T^"^ are simply convenient labels and do not correspond exactly to the retarded or 
advanced 0- functions of time as, for example, in the case of the single particle propagator of Eq. H81|l . 



15 



One can now check that T^GqT^ = 0, when i ^ k: 

Gofl'Gof^Go = {Doi)U{Dod ® «) G,' {Dok)h{Dok) ® (dt) 
= {Doi)ii{njnkDl - rijUkD^) (g) fiidl Gq^ 

{mjDl - ninjDl)ik{Dok) O Ukdl = 0, (100) 

where we used Eq. (jS^ for Gq^, and the fact that 



[riiUj 



nkD- - ninjnkD1)ggg{ninjnkD''f, - niUjUkDl) = 0. (101) 



It is therefore clear that the solution to the Faddeev equation, Eq. (|7Hjl . has the form 
Xi = + Xf-, where Xj^ and Xf^ satisfy the independent equations 

XP = lf,^ + lf,^Y.^oX^^ (102a) 

k^i 

In a similar way we have that 

Gof^Gof^G^ = {DoMDo^) ® «) G,' {Dok)ik{Dok) ® (4) 
= {Doi)ii{njnkDl - rijnkDf) nid\ Gq^ 

{niUjDl - ninjDl)ik{DQk) ® Ukd^ 
= {Doi)ii{njnkDl) fiidl Go^{ninjDl)ik{Dok) ® fikdl 

and 



Gof^^Gof^Go = {Doi)U{Doi) ® «) G^' {Dok)ik{Dok) ® {dt) 
= {Doi)ii{njnkDl - UjUkD^) ® Uid^ Gq^ 

{niUjDl - ninjDl)ik{Dok) ® rikdl 
= {Doi)ii{-njnkDf) ^ fiidl GQ\-ninjDl)ik{Dok) ® Ukdl 
= (DoMDW ® «) {D^ok)ik{Dok) ® (dt). (104) 

Thus, apart from external two-body propagators (-Dot) and {Dok), all other propagators 
explicitly shown on the RHS of Eqs. (jlOHj) and ()104|1 are truncated to their retarded or 
advanced parts. To take advantage of this simplification, we iterate Eqs. ()102|) once to 
obtain 



= ^^'^ + ^ E ^'''■(^0 + E ^oX;Go)f;« (105a) 
X,^ = + ^ E ^^"(^0 + E G'oX;Go)T,'^^ (105b) 



4 



16 



where 

Gof^'-G, = (27r)353(p^ _ p^) {D,,)U{Dl,) ® «), (106a) 

GoT^Go = (27r)353(p^ _ p^) {DlMDor) ® «), (106b) 

Goff'^G^ = {27,f6\v- - P.) (^o.)t.(/^o.) ® (di), (106c) 

Gof.^^Go = (27r)353(p^ _ p^) (Z},"Jt,(Z}o.) ® {d^), (106d) 

and where the amphtudes X[ and satisfy the Faddeev equations 

X: = lf[ + ^f[J2GoXl (107a) 

X^ = ^f,^ + ^ftJ2GoX'k, (107b) 

which are simpler than Eqs. ()102|) in that they utihze input t matrices whose adjoining 
propagators are either all retarded or all advanced: 

Gof[ Go = (27r)353(p^ _ p^) {D^,MDQ ® (108a) 

Go f t Go = (27r)353(p^ _ p^) {D^JU{D^J ® «). (108b) 

Note, however, that 

ii = Vi + ^Vi{Doi)ii, (109) 

so that the internal propagators, (Doi), used in constructing the physical two-body t matrices 
ti, retain both retarded and advanced parts. With the spectator 6 function removed, one 
can invert the Go's in Eqs. fllOHj) with the help of Eq. (jSEI), to obtain 

i gggn'ju'k {E - uj[^^) U{E - z){E- Uijk) UiUjUk ggg 



TUE) = — dz ", ' '-^^ , ' (110a) 

fiE) = — r dz ^^^^^'^'^^^ ~ ^^'jfe^ ^'^^ ujjjk) njUjUk ggg 

' 2tx {E - z - u'.f^- ie){E - z - ujjk- ie){z - Ui- it) 

where uj^^^ = i^i + uj'j + tu^, Uijk = i^i + oJj + cu^, uj'^^ = u'j + uj'^, and ujjk = ujj + ujk- In 
turn, this result shows that the Go in Eq. ()107a|) can be replaced by (Gg), and the Go in 
Eq. ()107bjl can be replaced by (Gg). Similarly, replacements can be made in appropriate 
Go's in Eq. (jl05|) and Eq. (|106|) . so that our final equations for the amplitude Xj are: 

X,=Xf + X,^ (111) 

= Iff + \Y.^f'- ( (GS) + 5^(GS)X;(GS) ) T;«, (112a) 

Xf = \ff + ^ E^^" ( <^o) + E(^o)^;(G'S) ) T,"^^, (112b) 
k^i \ j^k / 
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where the amphtudes X[ and satisfy the Faddeev equations 



Xl = lf[ + ^f[Y,{GdXl (113a) 
X:=^-ft + ^f,^Y.(G-o)X^k- (113b) 



2. Single degree of freedom equations 

The projection properties of the input two-body t matrices of Eqs. (IllUj) enable one to 
ehminate the doubled degrees of freedom from the in-matter Faddeev equations, Eq. ()113p . 
To see this, we first note that the amplitude Xj^ of Eqs. pi2ap is determined by the quantity 

(G5)X;(G5) = G'^r^lr^2r^3X;nln2n3G'^ (114) 

which has the doubled-degrees of freedom Faddeev amplitude XJ completely surrounded by 
projection operators n. The matrix structure of n is given in Eq. ()A29|) as 

n = ^{uj)^\uj) (115) 

where ^ {u) is a column vector defined in Eq. ()A10|) for bosons and in Eq. ()A20|) for fermions, 
and ^^{uj) is the corresponding row vector. Thus amplitude is expressible directly in 
terms of the single-degree of freedom amplitude 

X; = ^t(^^)^t(^^)^t(^^) x; '^(a;i)^(c^2)^(cu3). (116) 
Similarly, the matrix structure of projection operator n, given in Eq. ()A30|) as 

n = ±'r{uj)y\uj), (117) 

allows one to express amplitude X^ of Eq. ()112b|) directly in terms of the single-degree of 
freedom amplitude 

X- = i/\J^)r\J^)i/\J^) x; r(a7i)r(^2)r(^3). (ns) 

Moreover, it follows from Eqs. (jll3p that XJ" and X" themselves satisfy Faddeev equations 

^\ = \tI^\tIY.^''^1 (119a) 

XI = \ft ^\ftY.G^Xl (119b) 



where, with the spectator b function removed, 

T[{E) = — dz— ^^f-^^ — (120a) 



27r 7_oo {E - z - iUjf^ + ie){E - z - ujjk + ie) {z - Ui + ie) 



i (E- uj',,) ifiE - z)(E- tOi^k) 
Tm = ^ dzj^ ^ ^^--P (120b) 
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and 



tf = ^t(^)^t(^^)i~.^(^.)^(^^)^ (121a) 

if = rt(^)rtK) u 'r{u;,)rM. (121b) 

We note that in Eqs. (jll9p all doubled-degrees of freedom have been ehminated. In particu- 
lar, the two-body input to these equations is found from the convolution integrals, Eqs. (jl2(J|) . 
involving the single-degree of freedom t matrices tf and tf of Eqs. (jl2ip . 

Although both tf and tf are single-degree of freedom quantities, they themselves are 
constructed from a 16-component doubled-degree of freedom t matrix ti. Although it may 
not be difficult to solve the 16-component equation, Eq. p09|) . to obtain ii, it is useful to 
note that this equation can be recast into a four-component (2 x 2) equation as follows. 
Writing Eq. ^IW\i as 

ii = Vi + -Vi {ujUkDl - rijUkD^) ii, (122) 
straightforward use of Eqs. ()A29|) and ()A30|) allows one to write 

1 ^ . 

ti = Vi + -ViDiU, (123) 



where 



(124) 



(125) 



V. SUMMARY 

Using the real-time formalism, we have formulated equal-time three-body equations that 
describe three identical particles interacting via pair-wise interactions at finite temperature 
and density. Starting with the four-dimensional field-theoretic description of the 3—^3 
Green function, equal-time three-body equations were derived without resorting to any ap- 
proximations beyond that of the assumption of pair-wise interactions. 

Our resulting in-matter three-body equations, Eq. (fTHj) for the general case, and Eq. pil|) - 
Eq. ()113p for the case of instantaneaous potentials and free-like dressed propagators, have the 
familiar Faddeev form, although they differ from the usual zero-density Faddeev equations 
in that they involve doubled degrees of freedom (inherent in the real-time formalism), and 
they utilize one-body thermal Green functions (which have retarded and advanced parts, 
and depend on both temperature and chemical potential). At the same time, the form of our 
equations is similar to that of other formulations of the in-matter three-body problem, even 
though all other formulations have apparently been done either at zero temperature, or for 
non-zero temperatures, using the imaginary time formalism. However, what distinguishes 
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our approach in an essential way from all other derivations of the equal-time in-matter 
three-body problem, is that we have managed to avoid any approximations in the equal- 
time 3D reduction of the original 4D field-theoretic formulation. Moreover, our resulting 3D 
equations remain practical in that the equal-time two-body t matrix in three-body space, 
Tj, which determines the integral equation kernel, is given in terms of the 4D two-body t 
matrix by a simple convolution integral - Eq. (j72j) . 
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APPENDIX A: PROPAGATOR MATRIX STRUCTURE 



In the doubled degrees of freedom formalism, the matrix structure of propagators is 
discussed in detail, for example, in Ref. i^j. Here we give only a brief summary. 



1. Bosons 

In the doubled degrees of freedom formalism, the boson propagator is given by 2^ 

\ _ ■ ( cosh a; sinhu; \ ( dJ' \ / coshcj sinhu; \ ( ^^\ 

^ ' \ sinh oj cosh y I —d"' J \ sinh uj cosh uj j ^ ' 

where 

d'= , \ . , d^ = ^ -> ^=^-^ (A2) 

— UJ + le jr — ijj — le 2m 

and 



sinhcj = a//b(^)> coshcu = a/TT/^H, = _ ^ ■ (A3) 

Denoting 

rr / \ f coshcij sinhu; \ / 1 \ , . 

^-(")=(sinhc. coshcj' ^-=(o-lj' ^^4) 

one finds that Ub is not unitary but satisfies 

UB{oo)gBUl{u)=gB. (A5) 

It follows that 

[df(p)]-' = -^(p^-u;)gB. (A6) 
Comparison with Eq. (@)) provides explicit expressions for n and n: 

1 \ . . ^ / 



n = Ub{uj) ^ n = Ub{uj) q _ J (A7) 



20 



One can check the following important properties of n and n: 

n + n = qb (A8a) 
ngB^ = ngsn = (A8b) 
ngBn = n, ngBn = n. (A8c) 

For the case of bosons, these relations will define what is meant by the projection properties 
of the operators n and n. 

Eq. ()A7|) implies another convenient way of expressing n and n, namely, as products of 
column and row vectors: 

n = ^BM^iM, n = -rBH>B^M (A9) 

where 

/ f/fi \ / coshcj \ [ U?.\ ( sinhcj \ 

are the column vectors and '^J, 'J^, are the corresponding row vectors. It follows that 

'^i gs = 1, gB n = -1. (All) 

2. Fermions 

In the doubled degrees of freedom formalism, the fermion propagator is given by j2 

df{p) = ^^'^ ^i^^ U U COS.. - sin^ \ 
\^ — sm cos uj J \y a J ysmuj cos on J 

where cT, d", u are as in Eq. ()A2|) . and 

sin a; = V/fM, cos a; = ^1 - /^(u;), /f(c^) = ■ (A13) 

6 ~t~ J- 



Denoting 



[/i. = . , A14 

^ ' ' — sm cos uj ' ^ ^ 



it's clear that Up is unitary: 

?7F(u;)f/|(c<;) = uU^)Uf{uj) = 1. (A15) 

It follows that 

[df{p)]-' = -z{p'-u)(^l J). (A16) 
Comparison with Eq. provides explicit expressions for n and n: 

n = UAuj)(ll) f/|M, n = f/^M ( J ) (A17) 
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One can check the following projection properties of n and n: 

n + n = 1 (Al8a) 
nn^nn^O (Al8b) 
nn — n, nn — n. (Al8c) 

As in the boson case, we can write n and n in terms of column and row vectors: 

n^^F{oj)'^F{^), n = y^HrJlH (A19) 

where 

^^(^) =\jf] = \ . ' M^) =[jf] = \ ' (A20) 

are the column vectors and '^p^Yp^ are the corresponding row vectors. It follows that 

= 1 , = 1 . (A21) 

3. General 

By defining 



and 



c/=P-. ^=|!:« (bosoncase) 

1^ gli;' 1^ c/p 1^ ^Uf y Yf (fermion case) 

one can write, for both bosons and fermions, 

d^{v)^iU9{'i .° (A24) 



[d/(p)]-i = -i(pO-a;)^, (A25) 

UgU^ ^U^gU^g, (A26) 
= [/(jO)t/t, n = C/,(0;)c/t, (A27) 



n + n = g (A28a) 
ngn = ngn = (A28b) 
ngn = n, ngn = n. (A28c) 



n 
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Also, for both bosons and fermions, 

n = ^'^^ ^^g'^ = \, (A29) 

while 

n = ±rr^ gY = ±1, (A30) 

with the + sign being for fermions and — sign being for bosons. 
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